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ABSTRACT 

Context. Observations of hydrogen B aimer lines in solar flares remain an important source of information on flare processes in the 
chromosphere during the impulsive phase of flares. The intensity profiles of optically thick hydrogen lines are determined by the 
temperature, density, and ionisation structure of the flaring atmosphere, by the plasma velocities and by the velocity distribution of 
particles in the line formation regions. 

Aims. We investigate the role of non-thermal electrons in the formation regions of Ha, H/3, and Hy lines in order to unfold their 
influence on the formation of these lines. We concentrate on pulse-beam heating varying on a subsecond timescale. Furthermore, 
we theoretically explore possibility that a new diagnostic tool exists indicating the presence of non-thermal electrons in the flaring 
chromosphere based on observations of optical hydrogen lines. 

Methods. To model the evolution of the flaring atmosphere and the time- dependent hydrogen excitation and ionisation, we used a 
1-D radiative hydrodynamic code combined with a test-particle code that simulates the propagation, scattering, and thermalisation of 
a power-law electron beam in order to obtain the flare heating and the non-thermal collisional rates due to the interaction of the beam 
with the hydrogen atoms. To not bias the results by other effects, we calculate only short time evolutions of the flaring atmosphere 
and neglect the plasma velocities in the radiative transfer. 

Results. All calculated models have shown a time-correlated response of the modelled B aimer line intensities on a subsecond 
timescale, with a subsecond timelag behind the beam flux. Depending on the beam parameters, both line centres and wings can 
show pronounced intensity variations. The non-thermal collisional rates generally result in an increased emission from a secondary 
region formed in the chromosphere. 

Conclusions. Despite the clear influence of the non-thermal electron beams on the B aimer line intensity profiles, we were not able on 
the basis of our simulations to produce any unambiguous diagnostic of non-thermal electrons in the line-emitting region, which would 
be based on comparison of individual Balmer line intensity profiles. However, fast line intensity variations, well-correlated with the 
beam flux variations, represent an indirect indication of pulsating beams. 

Key words. Sumflares - Radiative transfer - Hydrodynamics - Line formation - Line:profiles 

1 . Introduction mal (gradual) phase, the heating leads to the evaporation of chro- 

mospheric and transition region plasma into the corona which is 
In the context of interareting flare loop hard X^ay footpoint gradual i y fille d with relati vely dense (up to ~ 10 10 cm" 3 ) and hot 
sources (Hoyngetal. 1981; Hudson & Farnik 2002), all con- / 10 7 K x flara n i asma dCzavknwska et al .111 999V The, radiation 



(~ 10 7 K) flare plasma (Czavkows ka et alJ [T999). The radiation 
temporary flare models (Sturrock 1968; Kovv & Pneuman 1976; from the flare region is now dom inatedby soft-X rays, EUV, and 
|Shibata||199ft |Turkmani et al.||2005|; |Fletcher & Hudson||2008D, again by emission in the optical spectral lines, 
regardless of their nature, assign a fundamental role during the We concentrate on modelling the formation of optically thick 
flare energy release, transport and deposition to the high-energy hydrogen spectral lines Ha> m and Hy in the early phases of 
non-thermal particle beams. In the impulsive phase of flares, the solar flares by the means of numerical radiative hydrodynamics 
beams formed by charged particles are also guided from the combined with a test particle approach to simulate the propa- 
acceleration site (wherever it is located) downwards along the gation> scattering and energy loss of an electron beam with the 
magnetic field lines into the transition region, chromosphere and power . law spec trum and prescribed time-dependent energy flux 
possibly photosphere. At lower atmospheric layers due to the propagating through the solar atmosphere and depositing its en- 
high density of local plasma, their kinetic energy is efficiently ergy int0 the solar plasma In Ms context we address three main 
dissipated by Coulomb collisions, the corresponding regions are questions' 
rapidly heated, and dramatic changes of temperatu re and ionisa- 
tion o ccur. This results in explosive evaporation dDoschek et alj 1. Does rapidly varying electron beam flux manifest itself in 
Il996h . The manifestations of the early flare processes can be the Balmer line intensities? 

obser ved in the microwaves, soft and ha rd-X rays, and optical 2. How do the non-thermal particles in the Balmer lines forma- 

lines (jTandberg-Hanssen & Emslie||1988|). Later on in the ther- t ion regions influence the line profiles and intensities? 

3. Can an unambiguous diagnostic method be developed that 

Send offprint requests to: Kasparova is applicable to observations of Balmer lines recognising the 
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presence of the non-thermal particles in the line formation 
regions? 

Due to the complexity of simultaneously treating non-LTE 
radiative transfer in deep layers of the solar atmosphere and the 
hydrodynamics (radiative hydrodynamics), only a few attempts 
have been made to model the optical emission of flares. First 



model s of pu lse-beam he ating were developed by Canfield et al 



1984) and Fisher et al (119850 . Recently, lAbbett & Hawlev 
1999 ) and I Alfred et al. (2005) studied emission in several lines 



and continua using complex radiative hydrodynamic simulations 
of electron beam heating on a time scale up to several tens of 
seconds. 

In this paper we concentrate on fast time variations on a 
subsecond time scale. In previous works on this topic, plasma 
dynamics was neglected. Simplified time-dependent non-LTE 
simulations of Ha line were then performed e.g. by using a 
prescribed time evolution of a flare atmosphere from i ndepen- 
dent hydrodynamic simulations of pulse-beam heating (Heinzel 
) or by solving approximate energy equation (Ding et all 



1991 



2001 



Both results showed significant Ha line response to 
pulse-beam heating on subsecond time scales. Here, we solve 
1-D radiative hydrodynamics of a solar atmosphere subjected to 
a subsecond electron beam heating and study emissions in the 
Ha, HJ3, and Hy lines. 

The paper is organised as follows: Section [2 describes the 
numerical code and models of beam heating. Results of simu- 
lations concerning flare atmosphere dynamics and B aimer line 
emission are presented in Section [3] There, we also analyse sev- 
eral proposed diagnostic methods for recognising the presence of 
electron beams in line formation regions. Section |4] summarises 
our results. 

2. Model 

The model covers three important classes of processes whose 
importance was identified in flares: 

1. Propagation of charged, high-energy particle beams with 
power-law spectra and time-dependent energy flux down- 
wards through the solar atmosphere and their gradual ther- 
malisation due to the Coulomb collisions with the am bient 
plasma in the solar atmosphere (Bai 1982; Emslie 19781) . 

2. The hydrodynamic response of low-/? solar plasma corre- 
sponding to the energy deposited by the beam. 

3. Time evolution of the ionisation structure and formation of 
optical emission in the chromosphere and photosphere where 
non-LTE conditions apply. 

The individual classes of flare processes are modelled using 
three computer codes, each modelling one class of the processes 
identified above. The codes have been integrated into one radia- 
tive hydrodynamic code. 



2.1. Flare heating 

The flare heating caused by an electron beam propagating 
from the top of the loop located in the corona (s = 9.5 x 
10 3 km corresponding to T = 1 MK) down wards is calcu - 
lated using a test-particle c ode (TPC) based on lKarlick v (1990); 
Karlicky & Henoux (1992). We assume an electron beam with 
a power- law electron flux spect rum [electrons cm -2 s" 1 per unit 
energy] (Nagai & Emslie 1984) 



F(E,t) = g(t)F(E) = g(t) (6-2)- 
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Fig. 1. Time modulations of the beam flux. Solid lines show g(t), 
dotted lines denote i g(t')df = E to t(t)/F max , see Eqs. ([]]) and 

0. 

Table 1. Model parameters 



Model 


E tot [erg cm" 2 ] 


g(t) 


^max [erg cm" 2 s" 1 ] 


6 


H_TP_D3 


9x 10 1U 


trapezoid 


4.5 x 10 1U 


3 


H_TP_D5 


9x 10 10 


trapezoid 


4.5 x 10 10 


5 


H_3T_D3 


9 x 10 10 


3 triangles 


6 x 10 10 


3 


H_3T_D5 


9 x 10 10 


3 triangles 


6 x 10 10 


5 


L_TP_D3 


1.5 x 10 iU 


trapezoid 


0.75 x 10 iU 


3 


L_TP_D5 


1.5 x 10 10 


trapezoid 


0.75 x 10 10 


5 


L_3T_D3 


1.5 x 10 10 


3 triangles 


1 x 10 10 


3 


L_3T_D5 


1.5 xlO 10 


3 triangles 


1 x 10 10 


5 



where S is the power-law index, g(t) e (0, 1) is a function de- 
scribing the time modulation of the beam flux, F m2LX is the max- 
imum energy flux, i.e. energy flux of electrons with E > Eq 
at g(t) = 1 . In order to model the electron spectrum by the test 
particles the beam electron energy is limited by a low and a high- 
energy cutoff, Eq = 20 keV and E\ - 150 keV, respectively. We 
present results for two types of time modulation g(t)\ three sym- 
metric triangular peaks and a trapezoidal modulation (see Fig.Q]) 
and two total deposited energies E tot (see Table [B 



f 

Jo 



g(t)dt, 



(2) 



where t\ is the duration of the energy deposit. The energy fluxes 
^max have been chosen in such a way that for both time mod- 
ulations the total deposited energy E tot is the same. The model 
parameters are specified in Table [T] their values are consistent 
with common beam characteristics derived from hard X-ray ob- 
servations. 
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The TPC simulates the propagation, scattering and energy 
loss of an electron beam with a specified energy flux and power- 
law index as it propagates through partly ionised hydrogen 
plasma in the solar atmosphere. The losses of the beam elec- 
tron kinetic energy caused by Coulomb collisions due to electron 
and neu tral (h ydrogen) components of solar plasma are given by 
Emslie(1978) 



A£ e < 



27T6> 4 



-A(x + s)nuv Q AtB 



2ne A 



A£en = —A (1 - X^H^eA^B , 



(3) 
(4) 



but also to distinguish between the beam energy deposited into 
the electron and hydrogen component of solar plasma and there- 
fore to calculate the non-thermal contribution to the transition 
rates in hydrogen atoms which is the crucial point for the present 
study. The test-particle approach used here naturally takes into 
account propagation effects of the beam and time evolution 
of ionisation structure of the atmosphere. This leads us to a 
more realistic d escription of beam energy lo sses as comp ared to 
the approach of lAbbett & Hawlevl (1 1999b or l Alfred et ail (120051) 
who used an analytic heating function corresponding to a sta- 
tionary solution of bea m propagation through the atmosphere 
dHawlev&Fisherill994l) . 



where E is the kinetic energy of the non-thermal electron, v e is 
the non-thermal velocity, nu = n p + n n is the number density of 
equivalent hydrogen atoms, n p and n n are the proton and neutral 
hydrogen number densities, x = n p /nu is the hydrogen ionisa- 
tion degree and s - 1 .4 x 10" 4 accounts for t he contribution from 
the m etals to the plasma electron density dHeinzel & Karlicky 
Il992l) . The metal contribution to electron densities is critical 
around the temperature minimum where the hydrogen is al- 
most neutral. We account for it in this approximate way. On 
the other hand, we neglect the helium contribution which can 
reach maximum 20% of total electron density in hig her alti- 
tudes. The Coulomb logarithms A and A' are given by lEmslid 
(1 19781) and At B is a constant TPC timestep which has to be cho- 
sen to satisfy the condition that the total beam electron energy 
loss AE = AEqq + AE Qn per a timestep is negligible relative to its 
kinetic energy, i.e. AE/E <?c 1. 

The scattering of the beam due to Coulomb collisions is 
taken into account using a Monte-Carlo method combined with 
the a naly tical expressions for the cumulative effects described by 
[Bail (Il982l) . The relation between the mean square of the beam 
electron deflection angle (0 2 ) and the corresponding energy loss 
AE (which holds if (6 2 ) «c 1 or equivalently if AE/E «c 1) is 
given by formula 



(0 2 ) = 



AE 



7l + 1 



(5) 



where y L = 1 / ^1 - v 2 lc 2 is the Lorentz factor. The new elec- 
tron pitch angle Oq + A# is given by 



cos(#o + AO) = cos #o cos 6 S + sin # sin 6 S cos . 



(6) 



where #o is the original pitch angle at the beginning of the time 
step, 6 S is given by equation ([5]) and using a 2-D Gaussian dis- 
tribution. The distribution of the azimuthal angle is uniform, 
0G(O,2;r). 

The TPC in principle follows the motion of statistically im- 
portant number of test particles representing clusters of electrons 
in the time varying atmosphere which responses through the 1- 
D HD code and the non-LTE radiative transfer code to the flare 
heating by TPC. The test particles with a time-dependent power- 
law spectra are generated in the corona at the loop-top and at 
each timestep the positions, energies and pitch angles of the par- 
ticle clusters are calculated. The macroscopic energy deposits 
into the electron £ ee and neutral hydrogen £ en component of so- 
lar plasma are obtained by summing the energy losses (AE QQ and 
AE Qn ) of a huge number of particle clusters for each position in 
the atmosphere using a fine equidistant grid. This approach al- 
lows not only to calculate the total flare heating 

*H = Sqq + S en 



2.2. 1 -D plasma dynamics 

The state and time evolution of originally hydrostatic low-/? 
plasma along magnetic field lines is calculated using a 1-D hy- 
drodynamic code. The temperature, density and ionisation pro- 
files of t he initial atmosphere correspond to the VAL 3C atmo- 
sphere (iVernazza et al.| [l981) with a hydrostatic extension into 
the corona. The half-length of the loop is 10 Mm. The time evo- 
lution of the atmosphere is initiated by the energy deposited by 
the beam. The main processes that determine plasma evolution 
in flare loops are: convection and conduction (both in 1-D due 
to the magnetic field), radiative losses and indeed the dominant 
factor is the flare heating here calculated by the TPC. The evolu- 
tion of plasma in the flare loop can be described by a system of 
hydrodynamic conservation laws 



o. 



dpu d , o dP ^ 
dE d d <9 

— + —{uE) = - T (uP) - -r c + S , 

Ot OS OS OS 



(7) 



(8) 



(9) 



where s is the position and u the macroscopic plasma velocity 
along the magnetic field line and p is the plasma density. The 
gas pressure and the total plasma energy are 



P = n n (l + x + e)k B T, 



E = 



P 1 2 
7 + , 

y - 1 2 



(10) 



where y = c p /c Y = 5/3 is the specific heats ratio, &b the 
Boltzmann constant. The hydrogen ionisation degrees x in the 
photosphere and chromosphere is calculated at each timestep by 
the time-dependent non-LTE radiative transfer code. In the tran- 
sition region and corona we assume x = 1 . The source terms on 
the right hand sides of the system of equations are: F g the paral- 
lel component of the gravity force in respect to the semicircular 
field line, Tc the heat flux, calculated using the Spitzer's classical 
formula, and 

S = <H-<R + Q 

includes all other considered energy sources and sinks, i.e. the 
dominant flare heating *K which drives the time evolution of the 
atmosphere, the quiet heating Q assuring the stability of the ini- 
tial quiescent unperturbed (hydrostatic) atmosphere and the ra- 
diative losses <R. Th e radiative losses are calculated according to 
Rosner etaD (1 19781) for optically thin regions and according to 
Peres et all (1 19821) for optically thick regions. 

The 1-D gas dynamics i s treated using the explicit LCPFCT 
solver (lOran & Borislll987l) . the Crank-Nicolson algorithm for 
the heat transfer and the time step splitting technique to couple 
the individ ual source terms of the energy equation with hydro- 
dynamics (lOran & Boris 20(X)|). 
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2.3. Time-dependent non-LTE radiative transfer 



3. Results of flare simulations 



Using the instant values of T, nu, and S en obtained by the hy- 
drodynamic and test-particle codes, a time-dependent non-LTE 
radiative transfer for hydrogen is solved in lower part of the loop 
in a 1-D plan-parallel approximation. The hydrogen atom is ap- 
proximated by a five level plus continuum atomic model. 

The level populations rii are determined by the solution of 
a time-dependent system of equations of statistical equilibrium 
(ESE) 

dm 



at 



(ID 



where contain sums of thermal collisional rates dj and radia- 
tive rate s Rg, and Rg are precondit ioned in the frame of MALI 
method (Rybicki & Hu mmed 1 19911) . The excitation and ionisa- 
tion of hydrogen by the non-thermal electrons from the beam 
are also included into Pij usin g the non-thermal collisional rates 
C nt following the approach of Fa ng et all (1 1993b 



C* = 1.73xl0 10 ^, 
=5.35xl0 9 ^, 



=2.94xl0 10 ^ 
ni 



C£ = 1.91x10*— . 



-^9 yen 



For transitions from the ground level we thus get 



Pij = Ri 



■ + C lj + C\). 



(12) 



(13) 



Non-thermal collisional rates from excited levels as well as 
three-body recombin atio n rates are not cons idered here since 
iKarlickvet all (12004 and lStepan et al.1 (l2007h found them to be 
negligible compared to thermal ones. 

In order not to bias the effects of the non-thermal collisional 
processes by effects caused by macroscopic plasma velocities, 
we excluded the advection term d(riiu)/ds from Eq. (fTTb . The 
omission of the advection term can be justified by small ve- 
locities (~ 10 km s" 1 ) in the Balmer line formation regions 
dNejezchleball 19981) attained during the first few seconds of the 
flare atmosphere evolution. On the other hand, the benefit is a 
significant simplification of the radiative transfer code. The sys- 
tem of ESE (fTTb is closed by charge and particle conservation 
equations 



(14) 



7=1 



where n Q is the electron density. The contribution of helium 
to ionisation is neglected. Because the electron density is not 
known in advance, the system of preconditioned ESE is non- 
linear due to products of atomic level populations (includ- 
ing protons) with n Q or n 2 Q . Therefore, the ESE and con- 
servation conditions (ITU) are linearised with respect to the 
level populations and electron density. The complete system 
of equations is then solved using the Crank- Nicolson algo- 
rithm and Newton-R aphson iterative method (Heinzel 1995; 
iKasparova eTai1l2003l) . 

The non-LTE transfer is solved on the shortest time step 
given by the time step splitting technique in the hydrodynamic 
part, see Section 12.21 Resulting electron density (ionisation) is 
then fed back to the hydrodynamic equations and the TPC. 



We computed the atmosphere dynamics and time evolution of 
the Ha, H/3, and Hy line profiles resulting from a time-dependent 
electron beam heating of an initially hydrostatic VAL C atmo- 
sphere. 



3.1. Flare dynamics 

Figures [2] and [3] show the time evolution of temperature, hydro- 
gen ionisation, energy deposit 9i, and energy deposit to neutral 
hydrogen £ en of models specified in Table [T] Shortly after the 
beam injection at the loop top at t = s, the chromosphere is at 
t ~ 0.25 s heated mainly at heights between ~ 1000 - 2000 km. 
The temperature rise is higher for models with triangular peak 
modulation (3T) than for the trapezoidal one (TP) due to the 
higher energy flux injected into the atmosphere (compare the 
time evolution of total injected energy in Fig.Q]). Since the low- 
energy electrons are stopped higher in the atmosphere, the total 
energy deposit for steeper beams (higher S) is larger than for flat- 
ter beams at these heights (compare l~i at t = 0.25 s in Figs. [2] 
and [3]) and the temperature rise is most significant for models 
L_3T_D5 and H_3T_D5, see the first panel in Figs. [2 and [3l re- 
spectively. This is generally true for all times during the heating; 
the temperature above s ~ 1000 km is larger for larger 6 when 
comparing models with the same time modulation git) and F max . 
On the contrary, at lower heights s < 1000 km due to the larger 
heating for flatter beams, see Figs. [2 and [3 the temperature at 
those atmospheric layers rises more for models with S = 3 than 
with S = 5. 

In the low-flux models the heating leads to a gradual increase 
of temperature above ~ 1000 km. The steep rise of temperature 
from chromospheric to coronal values is shifted by about 200 km 
from the preflare height to s ~ 1900 km, see Fig. O Heating by 
a higher flux results in a secondary region of a steep temperature 
rise at s ~ 1400 km which is formed at t ~ Is, see Fig. [3] Due 
to the locally efficient radiative losses, the temperature structure 
in all models up to s ~ 2000 km follows the time modulation of 
the beam flux; i.e. it rises and drops - compare the temperature 
structure e.g. at t = 1.5 s and t = 2 s for 3T models and TP 
models in Figs. [2] and [3] and the time evolution of temperature at 
two selected heights in Fig.|U 

Similarly to the temperature evolution, ionisation increases 
above s ~ 1000 km shortly after the beam injection, see the mid- 
dle panel in Figs. [2] and [3 Again, due to beam flux modulation 
and the dependence of the energy deposit on S, the increase of 
ionisation at those layers is most significant for 3T models with 
S = 5. As the heating continues, the layers above ~ 1200 km 
become completely ionised (at t ~ 2 s and t ~ 1 s for low and 
high flux models, respectively - see thick lines in Figs. [2] and [3]) 
and ionisation does not change during further heating - see also 
Fig. IH 

On the contrary, lower heights, below s ~ 1000 km, exhibit 
most significant increase of ionisation for models with S - 3. 
The ionisation at these layers reacts to the beam flux time modu- 
lation more for the flatter beams - see Fig. [4] which also demon- 
strates that the relaxation of ionisation to preheating values lags 
behind the time evolution of temperature. 

This effect previousl y shown by iHeinze 3 ([1991) and 
Heinzel & Karlickv (1992) is due to time evolution of the ratio 
of the number of recombinations to photoionisations. Detailed 
behaviour of phot oionisations followed by photorecombinations 
was discussed bv lAbbett & Hawlevl {l999, Section 4). 
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500 1000 1500 2000 500 1000 1500 2000 500 1000 1500 2000 

s [km] s [km] s [km] 



Fig. 2. Temperature, ionisation, and energy deposit corresponding to the low-flux models (L) and their evolution in time (from 
top to bottom). Different line styles denote four low-flux models. Thin solid line shows initial VAL C temperature and ionisation 
structure. Left: Temperature. Middle: Ionisation. Thick lines denote the models with C nt , thin lines without C nt . Right: Total energy 
deposit *H (thick lines) and energy deposit to hydrogen £ en (thin lines). Only deposits corresponding to the triangular modulation 
models (3T) are displayed. 
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500 1000 1500 2000 500 
s [km] 



1000 1500 
s [km] 



2000 500 



1000 1500 
s [km] 



2000 



Fig. 3. Temperature, ionisation, and energy deposit corresponding to the high-flux models (H) and their evolution in time (from top 
to bottom). The notation is the same as in Fig. [2] 
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L_3T_D3 
L_3T_D5 



s = 1500 km 



s = 1000 km 



t [s] 



I 


I I I 

H 3T D3 
H 3T D5 




*'\/' '\ s = 1500 km 


I 


s = 1000 km - 
i I I I 



s = 1500 km 



. s = 1 000 km 



t [s] 



t [s] 



Fig. 4. Time evolution of temperature (left) and ionisation (right) at two heights s = 1000 and 1500 km. Top: low-flux models 
L_3T_D3 (dotted) and L_3T_D5 (dash-dotted). Bottom: high-flux models H_3T_D3 (dotted) and H_3T_D5 (dash-dotted). C nt were 
included. 



3.2. Influence of non-thermal collisions! rates 

To evaluate the influence of the non-thermal collisional rates, 
two separate runs with and without C nt (Eq. [12]) were made for 
each model in Table [TJ 

Taking into account C nt leads only to marginal changes of 
temperature and density structure of the atmosphere (up to 17%). 
On the contrary, hydrogen ionisation and emission in B aimer 
lines may significantly differ in models with and without C nt . 
Generally, the effect of C nt is stronger for models with larger 
^max or lower S. In the low-flux models, C nt significantly modify 
the time evolution of ionisation structure. They lead to a faster 
complete ionisation of layers above s ~ 1200 km and cause an 
increase of ionisation in the layers below, compare thin and thick 
lines in Fig.0 The influence of C nt in the high-flux models is lo- 
calised mainly in the layers below s ~ 1000 km where the flatter 
beams increase the ionisation. The upper parts of the atmosphere 
are affected by C nt only temporarily, till t ~ 0.5 s, when they 
contribute to the fast ionisation of those layers - see Fig.[3j 

Since C nt are directly proportional to the energy deposit on 
hydrogen 8 en (Eq. [J2]), their influence is strongly linked to the 
8 en as a function of height. Consequently, C nt affect B aimer line 
intensities according to their formation heights. That can be un- 
derstood in terms of the contribution function CF to the outgoing 



intensity 1% 



I CF A ds = I Tj A (s)e 



T A (S) 



ds. 



(15) 



where r\x is the emissivity and Tx is the optical depth. Figure [5] 
demonstrates the effect of C nt on the Ha line for a low and high- 
flux model. In high-flux models, C nt affect mainly the line wings. 
A new wing formation region appears at heights of maximum of 
energy deposit on hydrogen, lower 6 producing a stronger con- 
tribution to CF. In the low-flux models, C nt influence also the 
line centre due to change of ionisation of upper layers where 
e.g. the Ha line centre is formed. The optical depth in the line 
centre at these heights is decreased and the Ha line centre forma- 
tion region is shifted deeper. Similarly as in the high-flux mod- 
els, a new region of wing formation region again appears at the 
height of the 8 en maximum, however the dominant part of the 
wing emission still comes from the photospheric layers, i.e. from 
s ~ 100 km - see Fig.[5j 

3.3. Time variation of line intensities 

The intensity variations depend on the maximum beam flux. The 
low beam flux results in gradual increase of intensities (models 
L), whereas high beam flux (models H) causes rapid heating of 
the atmosphere and hence fast and larger increase of line intensi- 
ties - see Fig. [6] which compares time variations of Ha, Hj3, and 



Kasparova, J. et al.: Response of optical hydrogen lines to beam heating: I. Electron beams 

E ^ erq cm_3 if ,] io» 10- 10- ED fo° [erg cm "io ?_,] 1* 




Fig. 5. Ha contribution functions for high (top) and low (bottom) flux models for trapezoidal beam flux modulation. Black solid 
curves indicate the total energy deposit, white dashed curves denote the energy deposit on hydrogen. Gray lines display contours of 
optical depth r = 0.1, 1, 10. Top left: Model H_TP_D3 with C nt . Top centre: Model H_TP_D5 with C nt . Top right: Model H_TP_D3 
without C nt . Bottom left: Model L_TP_D3 with C nt . Bottom centre: Model L_TP_D5 with C nt . Bottom right: Model L_TP_D3 without 
C nt . Colour scale denotes log CF. All plots correspond to t = 1.8 s- maxima of line intensities, see Fig. [6] 



Hy for models L_3T_D3, L_3T_D5, H_3T_D3, and H_3T_D5. 
(Owing to use of the five level plus continuum atomic model, re- 
sults concerning Hy line should be regarded only qualitatively.) 
Since larger 6 results in larger heating of the upper parts of the at- 
mosphere, it leads to higher line centre intensities of Ha and H/3 
which are formed in the upper parts of the atmosphere. On the 
other hand, the whole Hy line and Ha and Hj3 wings are formed 
in deeper layers, therefore they show more prominent intensity 
variations for flatter beams (6 = 3)- compare D3 and D5 models 
in Fig.[6l 

Figure [7] demonstrates the time variation and changes in- 
duced by C nt of formation regions of the Ha line centre 
(A A = A), line wing AA = 1 A, Hy line centre (A A = A), and 
Hy line wing (AA = 0.6 A) for L_3T_D3 model. Shortly after 
the first beam injection into the VAL C atmosphere, at t ~ 0.1 s, 
line intensities decrease by a factor of ~ 2 and more - see Fig. [6] 
Such a dip appears due to a temporal increase in optical depth, 
r, which is caused by C nt . A similar decrease appearing later, at 
t ~ 0.4 s is present also in the models without C nt - compare CF 
with and without C nt in Fig. [7] That dip is a result of relative im- 



portance of t he thermal collis ional rates C\j. This behaviour was 
explained by Heinze 3 (119911) for the case of a 3 -level model of 
hydrogen as a consequence of steeper rise of second level pop- 
ulation, ?72, in comparison to n^ with time. Figure [7] also shows 
that a region of secondary Ha wing emission is formed shortly 
after the beam injection and its intensity slowly increases in time. 

In the case of the high-flux model, the emission from the sec- 
ondary formation region dominates and photospheric contribu- 
tion to the wing intensities diminishes; this behaviour is typical 
for all studied B aimer lines. As regards Hy line, formation of the 
Hy line centre is due to C nt completely moved from the photo- 
sphere to the layers above s ~ 1200 km. Contrary to the model 
without C nt , no emission from the photosphere contributes to the 
outgoing intensity. Furthermore, C nt again lead to the prominent 
secondary wing formation region at s ~ 1000 km which occurs 
in the model without C nt at much later time t ~ 2 s - see last 
panel in Fig. [7] 

Intensities of all three studied lines show a good correlation 
with the beam flux on a time scale of the beam flux variation, i.e. 
on a subsecond time scale. Depending on the amount of heating, 
time variations are caused by the time-dependent temperature 
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structure, e.g. Ha line centre in the high-flux models, and in- 
fluence of C nt . Line intensities which are affected by C nt (line 
centres in the low-flux models and line wings in the high-flux 
models) show more significant time variations. Maxima of line 
intensities lag behind the beam flux maxima, the time lags are 
generally larger for line wing intensities than for line centre in- 
tensities. Such a time lag is not a beam propagation effect but it 
is related to different time variations of electron density at dif- 
ferent heights. We are aware of the fact that velocities would 
cause asymmetry of lines and modify the line intensities. Here 
we concentrate on the beam influence on line formation. For de- 
tailed comparison with observed Balmer lines, velocity should 
be considered in calculating the line intensities. 

3.4. Diagnostic tools 

Having obtained the time evolution of Balmer line profiles for 
various electron beam parameters, we can search for observ- 
able signatures that could provide a method for diagnostic of the 
electron beam presence in the Balmer line formation regions. In 
order to distinguish between the flare energy transport by the 
non-thermal electrons and other agents (e. g. Alfven waves - 
Fletcher & Hudson (2008)), we consider a method suitable for 
such a diagnostics only if C nt lead to significant and systematic 
differences in measured quantities. 

3.4.1 . Intensity ratios 

Recently, Kasha pova et al.l (I2008I) reported on so-called side- 
lobes in Ha/H/5 intensity ratio (i.e. increased value of Ha/Uj3 at 
AA ~ 0.5 A with respect to other wavelength positions) observed 
in flare kernels associated with radio and hard X-ray bursts. They 
attributed the appearance of such sidelobes to the effects of non- 
thermal electrons. 

Using our simulations we are able to check whether the ob- 
served sidelobes are a feature related to the electron beams. An 
example of a sidelobe is shown in Fig. [8] (left). It appears as a 
local maximum at AA ~ 0.6 A and varies on a timescale similar 
to the beam flux modulation. Due to velocities, such sidelobes 
could be asymmetric with respect to the line centre. On the other 
hand, there are observations of almost symmetric Balmer lines 
and Ha/Uj3 r atios in the flare k ernels associated with hard X- 
ray emission ( Kotrc et al. 2008), thus our models which do not 
consider velocities in radiative transfer can be applied to such 
situations. 

Figure [9] shows the time evolution of Ha/H/3 at several AA, 
R a p(AA, t) (indicated also in Fig. [8]), for all considered models 
of the beam heating (see Table [T]). In these plots, a sidelobe at a 
certain A A would appear as an increased line above the others, 
varying according to the beam time modulation. 

To consider a sidelobe being caused by the electron beam, 
the sidelobe should correspond to a simulation with C nt included. 
For some beam parameters R a p(0.5 A, t) exhibits the reported 
sidelobe behaviour with the exception of the high-flux models 
with S = 3 - see Fig. [8] (right) or the corresponding third panel 
in the top and bottom row in Fig. [9]- where R a p(0.5 A, t) rapidly 
drops below R a p(AA > 0.5 A, t) and no sidelobe exists in R a p 
on a beam timescale. Ha/Hy ratio shows a similar behaviour 
to R a p whereas H/3/Hy ratio is much weakly sensitive to beam 
parameters. 

On the other hand, neglecting the effect of C nt , i.e. assuming 
other agents than electron beams for the flare energy transport, 
sidelobes at AA = 0.5 A are present at all models. Thus, the 



R = Ha / HB noC nt R = Ha / HB C nt 




Fig. 8. Intensity ratio R a p(AA, t) as a 2D function of wavelength 
and time for the H_TP_D3 model. Left: without C nt , right: with 
C nt . The black dashed line denotes time modulation of the beam 
flux. Thick solid lines represent R a p at selected A/1 = 0,0.5,1 A 
which are shown for all models in Fig. |9l For display purposes, 
the time evolution is shown from t = 0.8 s. 
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Fig. 10. Time evolution of Ha (solid), UJ3 (dashed), and Hy 
(dot-dashed) wavelength-integrated intensities. The modulation 
of the beam flux is shown as colour-changing broken line, the 
colour indicates the flux value. Thick lines denote models with 
C nt and thin without C nt . Total deposited energy E tot is the same 
for both models. 



observed sidelobes cannot be considered as a unique signature of 
beam presence in the atmosphere but they are probably related 
to an impulsive heating. 

Note that the maximal sidelobes from our simulations may 
appear, for some models, at wavelength positions slightly differ- 
ent from AA ~ 0.5 A. 

Furthermore, other kinds of intensity ratios, e.g. a rela- 
tive line intensity with respect to the line centre value such as 
R a (AA, t) = I(AA, t)/I(0 A, t), do not show any systematic differ- 
ence between the models with and without C nt either. Therefore, 
the intensity ratios do not provide a reliable diagnostic tool suit- 
able for analysing the presence of the electron beams in the 
Balmer line formation regions. 

3.4.2. Wavelength-integrated intensity 

Wavelength-integrated intensity Z to t (proport i onal t o equivalent 
width) was recently proposed by lCheng et al.l (120061) as a tool for 
diagnostics of the non-thermal effects in the solar flares. On the 
basis of static semiempirical models, they propose to judge the 




Fig. 6. Time evolution of the Ha, UJ3, and Hy line intensities in three selected wavelengths for each line. Solid lines denote models 
with C nt , dashed without C nt . The black solid line shows the beam flux time modulation in relative units. From top to bottom: 
L_3T_D3, L_3T_D5, H_3T_D5, and H_3T_D5 model. 
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Fig. 7. Time evolution of CF of Ua (AA = A), Ua (AA = 1 A), Hy (A^ = A), and Hy (AA = 0.6 A) for L_3T_D3 model. Gray 
lines denote contours of = (0.1, 1, 10). Top: with C nt , Bottom: without C nt . Colour scale denotes log CF. 
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Fig. 9. The intensity ratios R a p(AA, t) for three selected wavelengths AA. Solid lines denote models with C nt and dashed without 
C nt . The black solid line shows the time modulation of the beam flux. Top: trapezoidal time modulation. Bottom: modulation by 
triangular peaks. 



relative importance of thermal and non-thermal heating in flares detail the wavelength-integrated intensity 
by analysis of / to t(Har) and / to t(CaII 8542A) which show signifi- 
cantly different sensitivity to C nt . Using this idea, we analyse in +A/W 

/ tot = k \ [I (A, t) - I(A, t = 0)] dA (16) 
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of the studied B aimer lines for more general radiative hydrody- 
namic models, k is a normalisation to scale intensities, here we 
use line centre intensity at t - s. 

From our simulations it follows that maximum reached val- 
ues of 7 t ot are predominantly given by the total deposited en- 
ergy E tot (see Eq.[2]and Fig.[lO|). Detailed time evolution of the 
beam flux is reflected as "loopy structures" in / tot - 7 tot plots - 
see the top row of Fig.Qj]- and the time variation of 7 tot corre- 
lates with the beam flux variation, see Fig.[lOl However, there is 
no unique dependence of B aimer line 7 tot on F max , see the right 
panel of Fig. [TOl which shows gradual increase of 7 tot with local 
minima and maxima corresponding to the time modulation of 
the beam flux. Moreover, in the case of the high-flux models, 7 tot 
depend also on 6. As a consequence of increased wing emission 
for lower S, see Fig.[6j 7 tot of all studied B aimer lines reach large 
values for flatter electron spectra (compare centre and bottom 
panel in Fig.fTTT). 

Taking C nt into account leads to a significant increase of 
7 to t - see Figs. [TO] and [TU but similar increase can be caused 
by stronger heating by other mechanisms than electron beams. 
Due to this reasons, neither wavelength-integrated intensities of 
B aimer lines are good indicators of electron beam presence in 
the B aimer line formation regions. 

4. Conclusions 

Presented radiative hydrodynamic simulations revealed the com- 
plexity of the response of hydrogen B aimer lines to the electron 
beam heating. At the same time, they proved to be a very useful 
tool to obtain answers to questions raised in the Introduction. 

1. We showed that the B aimer line intensities do vary on beam 
flux variation time scales, i.e. on a subsecond time scale. The 
time variations are caused by time evolution of the temper- 
ature structure, electron density and influence of the non- 
thermal collisional rates. Depending on the amount of the 
beam flux, time evolution of line intensities may also exhibit 
both fast (pulse like), and gradual (e.g. an increase of in- 
tensity on a time scale larger than the beam flux time vari- 
ation) components, see e.g. the case of model L_3T_D3 in 
FigO S uch behaviour is k nown for the Ha line from obser- 
vations dTrottet et a l. 2000). Therefore, we conclude that the 
fast pulse-like variations seem to be a good indicator of the 
particle beams, namely when correlated with HXR or radio 
pulsations. 

2. Influence of the non-thermal rates on the B aimer lines de- 
pends on the beam parameters, both the energy flux and 
power-law index. C nt significantly alter the ionisation struc- 
ture, leading to a modification of the line formation regions 
which are not ionised due to the heating. Depending on the 
beam parameters, C nt can affect line centres, wings or both, 
but generally C nt result in an increased emission from a sec- 
ondary formation region in the chromosphere. 

3. Concerning the diagnostic tools based on Balmer lines, ex- 
cept for the close correlation of the time variation of the 
beam flux and the line intensities, we did not found any sys- 
tematic behaviour that would uniquely indicate the presence 
of the non-thermal electrons in the atmosphere solely from 
observations of Balmer lines. Complementary information 
such as hard X-ray emission or spectral lines having different 
sensitivity to C nt , e. g. Ca II (8542 A), are needed to assess 
the presence of the non-thermal particles. 

In this model study, we analysed the influence of the beam heat- 
ing and the non-thermal collisional rates on hydrogen Balmer 
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Fig. 11. Time evolution of 7 tot (Ha), 7 tot (H0), and 7 tot (Hy). The 
colour corresponds to the beam flux value - see also colouring 
in Fig.[l0l Arrows indicate the direction of time evolution. Solid 
lines denote models with C nt , dashed lines the models without 
C nt . Top: H_3T_D3 model. Centre: H_TP_D3 model. Bottom: 
H_TP_D5 model. Total deposited energy E tot is the same for all 
panels. 

lines in the case of prescribed fast beam flux modulation. The 
next step is to compare the observed line emission with the sim- 
ulated one using the non-LTE RHD models for beam parameters 
inferred from hard X-ray or radio emission. In this way, the role 
of different flare energy transport mechanisms e.g. such as al- 
ternative heatin g of the chromosphere by Alfven waves recently 
proposed by iFletcher & Hudson] (120081) can be adequately ad- 
dressed. We plan to apply our code to fast time variations of 
Ha and hard X-ray emis sions observed during solar flares (e.g. 
iRadziszewski et al.ll2007l) . 
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